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Abstract 

As toy models for space-time on the Planck scale, we consider examples of fermion 
systems in discrete space-time which are composed of one or two particles defined on 
two up to nine space-time points. We study the self-organization of the particles as 
described by a variational principle both analytically and numerically. We find an 
efi'ect of spontaneous symmetry breaking which leads to the emergence of a discrete 
causal structure. 



Contents 

1 Introduction Q 

2 Preliminaries and Basic Definitions 

3 One-Particle Systems 

4 Two-Particle Systems in the Critical Case |^ 

4.1 The Local Correlation Matrices and the Fermion Matrix 11 

4.2 Two Space-Time Points 15 

4.3 Three Space-Time Points 17 

4.4 Four Space-Time Points: Spontaneous Breaking of the Parity Symmetry . . 21 

4.5 Five Space-Time Points: Spontaneous Breaking of the Translation Symmetry 24 

4.6 Six to Nine Space-Time Points: Emergence of a Two-Dimensional Lattice . |26 



Two-Particle Systems for the Variational Principle with Constraint 

5.1 Two Space-Time Points 

5.2 Three Space-Time Points: Emergence of a Discrete Causal Structure . . 



27 
2S 



A Numerical Implementation of the Critical Case |31 



B Numerical Implementation of the Variational Principle with Constraint [33 



References 36 



1 



1 Introduction 



Following ideas from [3], fermion systems in discrete space-time were introduced in [2] as 
systems of a finite number of quantum mechanical particles defined on a finite number of 
space-time points (see also the review articles O [6] ) . The interaction of the particles is 
described by a variational principle, where we minimize an action defined on the ensemble 
of the corresponding wave functions. A-priori, there are no relations between the space- 
time points; in particular, there is no nearest-neighbor relation and no notion of causality. 
The idea is that these additional structures should be generated spontaneously. More 
precisely, in order to minimize the action, the wave functions form specific configurations; 
this can be visualized as a "self-organization" of the particles. As a consequence of this 
self-organization, the wave functions induce non-trivial relations between the space-time 
points. The hope is that these relations give rise to additional structures, which, in a 
suitable limit where the number of particles and space-time points tends to infinity, can 
be identified with the local and causal structure of Minkowski space or a Lorentzian 
manifold. In this limit, the configuration of the wave functions should go over to a Dirac 
sea structure. 

This intuitive picture has been made precise to some extent in the following papers. 
In [4] the variational principle in discrete space-time is analyzed, and it is shown that min- 
imizers exist. In [7] the symmetries of discrete fermion systems are studied systematically. 
It is proved under general assumptions, which include all cases of physical interest, that 
the permutation symmetry of the space-time points cannot be respected by the wave func- 
tions. In other words, the permutation symmetry of discrete space-time is spontaneously 
broken by the particles. This shows that the wave functions indeed induce non-trivial 
relations between the space-time points. The paper [8] is devoted to the transition from 
discrete space-time to Minkowski space. The method is to analyze regularizations of vac- 
uum Dirac seas in composite expressions which arise in our variational principle. Finally, 
in [3j we study the so-called continuum limit. Simply assuming that the wave functions 
form a Dirac sea configuration, we get an effective continuum theory where Dirac particles 
interact with classical gauge fields. Our variational principle determines the corresponding 
gauge groups and the coupling of the gauge fields to the Dirac particles. We thus obtain 
many promising concrete results which seem to indicate that our variational principle 
really is of physical significance. 

The motivation for the present paper is that the representation theory in [7j does not 
give an intuitive understanding of how the mechanism of spontaneous symmetry breaking 
works. Also, the abstract constructions do not give information on what the resulting 
smaller symmetry group is. More specifically, we would like to know whether we get 
relations between the space-time points which can be interpreted as a discrete version 
of a causal structure. In order to get a better understanding of these important issues, 
we here study concrete examples involving one or two particles. Clearly, these systems 
are much too small for being of direct physical interest. But they can be considered as 
interesting toy models, where the general mechanism of spontaneous symmetry breaking 
and the generation of a discrete causal structure can be studied in detail. 

This paper begins with the basics and can be used as an easily accessible introduction to 
fermion systems in discrete space-time. Whenever possible, we give explicit calculations 
and self-consistent proofs. But we also present numerical results and refer to general 
theorems in jUI?], which instead of reproving we illustrate or explain in words. 
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2 Preliminaries and Basic Definitions 



In this section we give the basic definitions and cohect a few general properties of fermion 
systems in discrete space-time. For simphcity, we always restrict attention to the case of 
spin dimension one, because all our examples will be constructed in this case. For the 
generalization to higher spin dimension we refer to [HIT] • Let H he a finite-dimensional 
complex vector space, endowed with a sesquilinear form <.|.> : H x H ^ C, i.e. for all 
u,v,w G H and a, (3 £ C, 

<u I av + f3w> = a <u \ v> + f3 <u \ w> 
<au + f3v\w> = a <u\ w> + (3 <v \ w> . 

We assume that <.|.> is Hermitian symmetric, 

<u I v> = <v I u> , 

and non-degenerate, 

<n I i;> = \/veH =^ u = 0. 

We also refer to {H, <.|.>) as an indefinite inner product space. We point out that the inner 
product does not need to be positive, but it can have a general signature {p, q) (see [H |9] 
for details). Nevertheless, many constructions from Hilbert spaces can be carried over. In 
particular, we define the adjoint of a linear operator A : H ^ H hy the relation 

<u\Av> = <A*u\v> yu,veH. 

A linear operator A is said to be unitary if A* = A^^ and self-adjoint if A* = A. It is 
called a projector if it is self-adjoint and idempotent. 

Let M be a finite set consisting of m = points. To every point x £ M we associate 
a projector E^- We assume that these projectors are orthogonal and complete in the sense 
that 

ExEy = 5^yE^ and ^ = 1 . (2.1) 

Furthermore, we assume that the images Ex{H) <Z H oi these projectors are non-de- 
generate subspaces of which all have the same dimension two and the same signa- 
ture (1, 1) (we remark that in the generalization to higher spin dimension, these subspaces 
have signature (n, n) with n > 1). The points x S M are called discrete space-time 
points, and the corresponding projectors Ex are the space-time projectors. The struc- 
ture {H,<.\.>.,{Ex)x&m) is called discrete space-time. 

To describe the particles of our system, we introduce one more projector P on H, 
the so-called fermionic projector, which has the additional property that its image P{H) 
is negative definite. The vectors in the image of P can be interpreted as the quantum 
mechanical states of the particles of our system, and thus we call the rank of P the 
number of particles, and denote this number by / = d\m.P{H). The name "fermionic 
projector" is motivated from the fact that in physical applications [3], the particles are 
Dirac particles and thus fermions. We refer to {H, <.|.>, {Ex)x<^m, P) as a fermion system 
in discrete space-time. 

A space-time projector Ex can be used to restrict an operator to the subspace Ex{H). 
We refer to this restriction as the localization at the space-time point x. Localizing the 
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fermionic projector gives rise to the so-called discrete kernel P{x,y) := E^PEy, which 
maps the subspace Ey{H) C H to Ex{H) and vanishes otherwise. It is most convenient 
to regard it as a mapping only between these subspaces, 

P{x,y) : Ey{H) ^ Ex{H) . 

Then the product of P{x,y) and P{y,x) is an endomorphism of Ex{H), 

A,y = Pix,y)P{y,x) : E,{H) ^ E,{H) , 

referred to as a closed chain. The closed chain is a self-adjoint operator on Ex{H). The 
roots of its characteristic polynomial counted with multiplicities are denoted by A+ and A_. 
We define the spectral weight of Axy by 

\-^xy\ = |A+| -|- |A_|. 

Similarly, the spectral weight of the square of the closed chain is defined by \A^y\ = |A+p-|- 
|A_p. Summing the spectral weights over all space-time points we get a positive number, 
which only depends on the fermionic projector relative to the space-time projectors. For 
a given k > we consider the family of fermionic projectors 

V{k) := \ P with ^ \Axy\^ = k > (2.2) 
[ x,y&M J 

and introduce the function 

Z[P] = \Aly\ . (2.3) 

x,y&M 

Our variational principle is to 

minimize Z[P] by varying P in V{k) , (2-4) 

keeping the number of particles / as well as discrete space-time fixed. For clarity, we 
sometimes refer to ()2.4p as the variational principle with constraint. 

Using the method of Lagrange multipliers, for every minimizer P there is a real pa- 
rameter /X such that P is a stationary point of the so-called action 

S^[P] = Yl ^M-y\ (2-5) 

with the Lagrangian 

^lA-^xy\ '■— \-^xy\ ~ l^a;?;! • i'^-^) 

A possible method for constructing stationary points is to minimize 5^, and this leads us 
to the so-called auxiliary variational principle 

minimize 5^[-P] by varying P . (2.7) 

Without particles (/ = 0), the variational principles (|2.4p and (|2.7p are of course 
trivial. If on the other hand there is only one space-time point (m = 1), there is only one 
closed chain A = P^ = P and hence 

l^l^l = 1^1 = /, = f-iif. 
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Thus the variational principles (j2.4p and (j2.7p are again trivial. Therefore, we shall always 
assume that 

/ > 1 and m > 1 . 



Before discussing the above variational principles, it is convenient to bring the fermion 
system into a more explicit form by writing down the projectors (-E'a;)agAf and P as matrices 
in a given basis of H. We shall always work with a basis of the following form: In each 
subspace Ex{H) C H we choose a pseudo-orthonormal basis (ef,e|), i.e. 

<ef \ej> = Si 5ij with si = 1 , S2 = — 1 . 

The orthogonality relation in (12. ip yields that <e^|ej> = for all x ^ y and all i,j G 
{1,2}. Furthermore, the completeness relation in (|2.ip yields that the vectors {Gf )^^l^2} 
form a basis of H. In particular, we see that H has dimension 2m and signature (m, m). 
Ordering the basis vectors as {e\, e\, e\, e\, . . .) , the inner product of H has the represen- 
tation 

<u\v> = {u\Sv) yu,v£H, (2.8) 

where the round brackets denote the standard Euclidean scalar product on C^™, and the 
so-called signature matrix S in block matrix notation is given by 




In the same block matrix notation, the space-time projectors become 




(2.10) 



It is instructive to briefly discuss the fermionic projector and the closed chain in this basis. 
Since P{H) is negative definite, the number of particles must be bounded by the negative 
signature of H, i.e. we always have f < m. The self-adjointness of P can be expressed by 

(P5)t = PS, (2.11) 

where the dagger denotes the transpose, complex conjugate matrix. For any x,y €z M, the 
closed chain Axy is a self-adjoint operator on Ex{H). Thus it is expressed by a (2 x 2)- 
matrix satisfying the relation 

Applying the identity det(i?C — Al) = det(C-B — Al) to the characteristic polynomial 
of Axy , we find that 

det(^^y-Al) = det(P(x,y)P(y,x) - Al) 

= det{P{y,x)P{x,y) - XI) = det{Ayx - XI) , 

proving that the operators A^y and Ay^ have the same characteristic polynomial. This 
shows that our Lagrangian is symmetric in the two space-time points, 

£^[Axy] = C^yx] yx,yeM. (2.12) 
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In contrast to Hilbert spaces, the self-adjointness of a closed chain A does not imply that 
its eigenvalues are real. But the calculation 

det{A - Al) = det(A^-x) = det( s{A^ -X)s) = det(^*-A) = det(A-Al) 



shows that its characteristic polynomial has real coefficients. This means that the \± are 
either both real, or else they form a complex conjugate pair. 

We now give a brief overview of the existence theory. In [3] it is proved that the 
variational principle with constraint (j2.4h is well-posed. 



Theorem 2.1 //'P(k) is non-empty, the variational principle 1^2. 4\) attains its minimum. 



The auxiliary variational principle behaves differently depending on the value of n, as the 
following consideration shows. 

Proposition 2.2 If fi < \, the Lagrangian L^j, is non-negative. If conversely fi > ^ 
and f >2, the action 5^ is not bounded from below. 

Proof. The Cauchy-Schwarz inequality 2|A+A_| < |A+p + |A_p yields 

|^|2 = (|A+| + |A_|)^ = |A+|2 + 2|A+A_| + |A_|2 < 2|A+|2 + 2|A_|2 = 2\A'^\ . 

Hence 

C, = \A'\-fi\A\' > Q-^ 

proving that the Lagrangian is non-negative in the case fJ. < ^. 

In the remaining case /x > ^, we shall construct a family of fermionic projectors 
where 5^ tends to minus infinity. In preparation, we briefly consider the simple case 
when H has signature (1, 1). In this case, the fermionic projector must be of rank one. It 
is most conveniently represented by choosing a vector u G H with <u | u> = — 1 and to 
set 

P = —luXul = —u^u^S, 



where we used a bra/ket-notation and then represented the inner product in the form (12. Sp 
with S = diag(l,— 1). Choosing u = {^/a, y/a^\^)'^ with a parameter a > 0, we obtain 
the fermionic projector with the matrix representation 

-a (3 



P = with/3 = Va(a + 1). (2.13) 

Coming back to the situation of the proposition, we know that / > 2 and thus m > 2. 
Working in the matrix representation (j2.9l 12.10]) . we choose a fermionic projector which 
is invariant on the subspace Ei{H) © E2{II), and on this subspace has the form 



P\EiiH)®E2{H) 



(-a 13 \ 

a + l -(3 

/3 -a 

\-l3 a + lj 



Since this matrix is built up of components of the form (j2.13p . it is obviously a projector 
on a two-dimensional, negative definite subspace. Thus, choosing P on the orthogonal 
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complement {Ei{H)®E2{H))-^ equal to any fixed projector on a negative definite subspace 
of dimension / — 2, we obtain a fermionic projector. 
Computing the closed chains, we obtain 

^11 = A22 = a^l + 0{a), Ai2 = A2i = -a^l + 0{a). 

whereas all other closed chains are independent of a. Hence 

\Aly\ = 4a^|lc2|+0(a3), ^ = Aa'' \\cA^ + 0{a^) . 

x,yGM x,yGM 

Using that \l£2\ = 2, we obtain for the action the asymptotic formula 

= 8a^(l-2/i) + 0(a='), 

and this really tends to minus infinity as a ^ cxd. ■ 

According to this proposition, in the case > ^ the auxiliary variational principle is 
not well-posed (except if there is only one particle, see Proposition 13.21 below). In the 
case /U < ^ a general existence theorem is proved in [4]. 

Theorem 2.3 If fi < the auxiliary variational principle attains its minimum. 

In the remaining so-called critical case /i = ^ of the auxiliary variational principle, partial 

existence results are given in [4], but the general existence problem is still open. From 

the physical point of view, the critical case is the most interesting case, because it allows 

to model a system of massive Dirac seas (see [El |3]). For this reason, we introduce the 

abbreviation C = Ci and S = Si. Using the transformation 
2 2 

\A'\-l\Af = (|A+p + |A_n-i(|A+| + |A_|f = ^(|A+|-|A_|f , 

the critical Lagrangian can be written in the simple form 

C[A] = i(|A+|-|A_|f . (2.14) 

Thus one can say that the critical case of the auxiliary variational principle tries to achieve 
that the X± have the same absolute value. 

We next discuss the symmetry structure of discrete fermion systems (for details see [7j ) • 
It is useful to distinguish between inner and outer symmetries. Inner symmetries are 
also called gauge symmetries and keep M fixed, whereas outer symmetries also involve 
a transformation of the space-time points. A unitary mapping U of is called a gauge 
transformation if it leaves the space-time projectors unchanged, i.e. 

UE^U'^ = E.J, Vx G M. 

The fermionic projector behaves under gauge transformations as follows, 

P — > UPU'^ . (2.15) 
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If P = UPU^^, the transformation U describes a gauge symmetry. The group of ah gauge 
transformations is denoted by ^. In the basis l\2.9\ I2.10p . a gauge transformation U can 
be written as 

/C/i \ 

U = •.. with [/, e U{l,l) . (2.16) 

VO uj 

This means that a gauge transformation sphts into the direct sum of unitary transforma- 
tions which act "locally" on the spaces Ex{H). We also see that Q is isomorphic to the 
group [/(l,!)'". 

To describe outer symmetries, we permute the space-time points and demand that, 
after a suitable unitary transformation, the discrete fermion system should be unchanged. 
We denote the symmetric group of M (i.e. the group of all permutations) by Sm- 

Definition 2.4 A subgroup O of the symmetric group Sm is called outer symmetry 

group of the discrete fermion system if for every a £ O there is a unitary transforma- 
tion U such that 

UPU-^ = P and UE.xU-^ = E^^^^ VxEM. (2.17) 
We next specialize j7i Theorem 9.1] to the case of spin dimension one. 

Theorem 2.5 Let (H, <.\.> , {Ex)x£M , P) be a fermion system in discrete space-time, 
where the subspaces Ex{H) C H have signature (1,1). Suppose that the number of parti- 
cles f and the number of space-time points m lie in the range 

1 < f < m - 1 . 

Then the system cannot have the outer symmetry group Sm- 

This result has the interpretation that the permutation symmetry of discrete space-time is 
spontaneously broken by the fermionic projector. It implies that the fermionic projector 
induces non-trivial relations between the space-time points. Unfortunately, the theorem 
does not give information on how these additional structures look like. One would hope 
for a structure which can be regarded as a discrete analogue of a causal structure. Such 
a structure can indeed be defined in general. Recall that after (|2.12[) we saw that the A-t 
are either both real or else they form a complex conjugate pair. This distinction can be 
used to introduce a notion of causality. 

Definition 2.6 (discrete causal structure) Two discrete space-time points x,y £ M 
are called timelike separated if the roots X± of the characteristic polynomial of A^y are 
both real. Conversely, they are said to be spacelike separated if the X± form a complex 
conjugate pair. 

This definition was first given in [6], and in [U [8] it is shown that for Dirac spinors in 
Minkowski space it gives back the usual notion of causality. The definition of the discrete 
causal structure can also be understood from the fact that it reflects the structure of the 
critical Lagrangian ()2.14p . Namely, if x and y are spacelike separated, the roots \± of the 
characteristic polynomial of A^y form a complex conjugate pair. Hence |A+| = |A_| = |A_|, 
and thus the Lagrangian p.l4p vanishes. Computing first variations of the Lagrangian, one 
sees that these also vanish, and thus P{x, y) does not enter the Euler-Lagrange equations. 
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This can be seen in analogy to the usual notion of causality in Minkowski space where 
points with spacelike separation cannot influence each other. 

The next simple proposition illustrates our definition of a discrete causal structure. 
The method is worked out in a more general context in [4, Section 4] and is used in the 
existence proof. 

Proposition 2.7 Every space-time point has timelike separation from itself. 

Proof. For any x £ M, we can use the idempotence and self-adjointness of our projectors 
to write the expectation value of P{x, x) as follows, 

<u I P{x, x) u> = <u\ ExPEx u> = <u\ ExP'^Ex u> = <PExU \ PEx u> . 

Since the image of P is negative definite, it follows that 

<u\P{x,x)u> < \fueH. 

Expressed in the basis (|2.91 12.10]) . we conclude that the matrix —sP{x,x) is positive semi- 
definite on the standard Euclidean C^. This implies that its determinant is non- negative. 
Using that the determinant is multiplicative and that det(— 1^2) = 1 and dets = —1, we 
conclude that 

detP{x,x) < 0. (2.18) 

The operator P(x,x) is self-adjoint on Ex{H) C H. Following the consideration af- 
ter (I2.12p . the roots n± of its characteristic polynomial either are real or they form a 
complex conjugate pair. In the latter case, the determinant of P{x, x) would be strictly 
positive, in contradiction to (I2.18|) . We conclude that the ^± are real. Since 

the spectral calculus yields that A-t = /x^ S M. ■ 

The aim of the present paper is to specify the spontaneous symmetry breaking of 
Theorem 12.51 and to clarify the connection to the discrete causal structure of Definition 12.61 
in simple examples. 

3 One-Particle Systems 

In this section we will show that for systems of only one particle, our variational principles 
can be solved analytically. We first remark that, since P has rank one, for any x,y €z M 
the discrete kernel P{x, y) has rank at most one. Thus the closed chains Axy also have 
rank at most one, and thus 

\Aly\ = \Axy\^ = TiiAxy)^ yx,yeM. (3.1) 

This implies that the functional in the variational principle with constraint (12. 4p is equal 
to K, making the variational principle trivial. 

Proposition 3.1 In the case f = 1 of one particle, every P G V{k) is a minimizer of the 
variational principle \2.4^ . 
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Using (|3.1|) . the Lagrangian of the auxihary variational principle (|2.6p simplifies to 

To further simplify the action, we introduce for any x G M the local trace px by 

Px = TtiExP) G M . (3.2) 

Next we choose a vector n 7^ in the image of P and normalize it to <u \ u> = — 1. Then 
the fermionic projector can be written in bra/ket-notation as 

P = -|u><n| . (3.3) 

Furthermore, the local trace becomes px = —<u\Ex\u>, and thus 

TT{Axy) = Tr{ExPEyP) = <u\Ex\u><u\Ey\u> = px Py ■ 

This simple factorization allows us to write the whole action as a square, 



S^ = Tr(^-y)' = (1-/^) E p'pI = (l-/")fE/'') 

x,y&M x,y£M \xeM / 

Now we can use the Cauchy-Schwarz inequality to obtain the following result. 



(3.4) 



Proposition 3.2 Consider the auxiliary variational principle \2. ?| ) for fermion systems in 
discrete space-time with one particle. Depending on the value of ji, there are the following 
cases: 

(i) If p = 1, Sf^ vanishes identically, and every fermionic projector is a minimizer. 

(ii) If p < 1, the minimum is attained. We can arrange by a gauge transformation that 

in the basis 112.91 \2.10\) . the minimizing fermionic projector can be written as ^3. 3\) 
with 

1 /O^ 



Exu = -j=[^^j VxGM. (3.5) 
(iii) If p > 1 and m > 1, the action is not bounded from below. 

Proof. The case (i) is obvious from (j3.4p . To prove (iii), we consider for a > the 
fermionic projector (j3.3p with 



^2n=(^^iL_^, EyU = Vy G {3,...,m}. 

A short calculation shows that 5^ as given by (13. 4p tends to minus infinity as a ^ 00. 
In the remaining case (ii), we apply the completeness relation in (j2.ip to (j3.2p . giving 

Y,Px = TrP = 1. 
Hence we can apply the Cauchy-Schwarz inequality to obtain 

1= (EpX ^ fEilfE^^l =-E/>^ (3-6) 

\xeM / \x£M / \x€M / x&M 
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Using this inequality in (|3.4p gives the lower bound 

1-M 



5„, > 



and equality holds if and only if the px are all equal. Hence every minimizer must satisfy 
the relations 

Px = — Vx G M . 

m 

This means that the components of the vector ExU =: (a, b)^ satisfy the relation a? — h'^ = 
— 1/m. Thus by an unitary transformation on Ex{H) we can arrange (j3.5p . ■ 

We remark that this proposition immediately generalizes to higher spin dimension. 

The minimizer in case (ii) has the property that the particle is completely delocalized 
in space-time, because according to ()3.5p it has the same probability 1/m to be at any of 
the space-time points. The resulting system is obviously permutation symmetric, and thus 
no spontaneous symmetry breaking occurs. Furthermore, as the closed chain Axy has rank 
at most one, one of the roots of its characteristic polynomial must vanish. Hence the A-t 
cannot form a complex conjugate pair. Following Definition 12. 6^ all pairs of space-time 
points are timelike separated, and the discrete causal structure is trivial. 



4 Two-Particle Systems in the Critical Case 

4.1 The Local Correlation Matrices and the Fermion Matrix 

We begin the discussion of many-particle systems with the simplest situation of two par- 
ticles (/ = 2) and the auxiliary variational principle in the critical case (j2.71 12.51 12.14p . 
For systems of two particles, the following construction is very useful for visualizing the 
discrete fermion system. The image of P is a two-dimensional, negative definite sub- 
space of H. Choosing an orthonormal basis (111,^2) (i.e. <Ui\uj> = —Sij), the fermionic 
projector can be written in bra/ket-notation as 

P = -|ni><ui| - \u2><U2\ . (4.1) 

For any space-time point x G M, we can introduce the so-called local correlation matrix 
Fx by 

Fx = {iFxyjkj=i,2 with (Fx)] = -<Ui\ExUj>. (4.2) 

The matrix Fx is Hermitian on the standard Euclidean C^. Thus we can decompose it in 
the form 

Fx = -{pxl + Vxff), (4.3) 

where a = (cr^, cj^, cr^) are the Pauli matrices, and px G M, Vx G M^. Taking the trace 
of (jOj) and using ([O]) . we find 

Px = - ^ <Ui\ExUi> = Tt{ExP) , 

i=l,2 

showing that px is again the local trace (|3.2p . The vector Vx describes the correlations of 
the two particles at the space-time point x. Using the terminology which is common in 
atomic physics and in quantum information theory for the description of q-bits, we refer 
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to Vx as the Block vector. Summing (j4.2p over x G M and using the completeness relation 
in (j2.ip . we find that X^xeA/ = ^ or, equivalently, 



Px 



and 



xGM 



x&M 



. 



(4.4) 



Furthermore, as the inner product in (j4.2p has signature (1,1), the matrix Fx can have 
at most one positive and at most one negative eigenvalue. Expressed in terms of the 
decomposition (|4.3p . this means that 



> Px 



Vx G M 



(4.5) 



We point out that, since in (|4.2p we took expectation values, the local trace and the 
Bloch vector are gauge invariant quantities. However, our construction clearly involves 
the arbitrariness of choosing the orthonormal basis (iii,U2) of P{H). More precisely, we 
have the freedom to unitarily transform the basis as follows, 



Ui 



with V G U{2) . 



(4.6) 



As is well-known from the transformation law for Pauli spinors in non-relativistic quantum 
mechanics (see for example |12l Section 3.2]), a unitary transformation of Pauli matrices 
corresponds to an S'0(3)-transformation of the vector index a G {1,2,3}, 



f3=l 



with 



R G S0{3) , (4.7) 
and thus the C/(2)-transformation (j4.6p leads to a joint rotation of all Bloch vectors, 

3 

< — VxGM. 

/3=1 

Hence we can say that the Bloch vectors Vx G are unique up to rotations (i.e. orientation- 
preserving isometrics) of the standard Euclidean M'^. 

For calculations, it is convenient to choose the basis (12.91 I2.10p and to introduce a 
matrix whose columns are the components of the orthonormal basis vectors ui and U2, 



/ (ni)i 
(^1)2 



(""2)1 \ 
(-"2)2 



(4.8) 



\{ui)2m {U2)2m) 



The matrix ^ is referred to as the fermion matrix. Clearly, the fermion matrix charac- 
terizes the fermion system completely. Gauge transformations (I2.16P act on ^ from the 
left, whereas the C/(2)-transformations (j4.6p act on it from the right. Thus we have the 
freedom to transform the fermion matrix according to 



with [/ G a and y G U{2) . 



(4.9) 
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Using the fermion matrix, we can write (j4.1|) in the compact form 

p = -^-^-t^- (4 10) 

For any x G M, the local correlation matrix becomes 

= -^^SE,^, (4.11) 

whereas the local trace and the Bloch vectors can be computed by 

= -Tii^^SE^^) , < = -TV(cj° ^^SE^^) . (4.12) 

Note that when applying the transformation (j4.9p in (|4.10p . the matrix V drops out, 
leaving us with the gauge freedom (j2.15|) . In (|4.11|) . on the other hand, the matrix U 
drops out, so that we are left with the C/(2)-freedom (j4.6p . 

We now collect a few results which clarify the significance of the local correlation 
matrices Fx- First, it is a useful fact that our variational principles can be expressed in 
terms of the local correlation matrices. 

Proposition 4.1 The roots X± of the characteristic polynomial of the closed chain A^y 
are written in terms of the local traces and the Bloch vectors by 



A± = ^ yPxPy + V^Vy ± ^J\pxVy + PyVx\^ - X UyPj . (4.13) 

Furthermore, 

A+A_ > 0. (4.14) 
Proof. Using (j4.10p . the closed chain can be written as the operator product 

Axy = Ex^^'^SEy^l^'^SEx. 

Using the identity det(i?C — Al) = det(C-B — Al), a cyclic commutation of the factors in 
the operator product does not change the spectrum. Hence the A-t are also the roots of 
the characteristic polynomial of the matrix 



^SEr^^^SE,,^ ^^S^ FrF, 



x-t y 



Using the decomposition (|4.3p of the local correlation matrix together with the iden- 
tities of Pauli matrices (see |12l Section 3.2]) 



3 
7=1 



(where e"^''' is the totally anti-symmetric Levi-Civita symbol), we obtain 

FxFy = ^ (PX + VxO) {Py + VyU) 

= \ {pxPy + VxVy + {pxVy + PyVx + iVx X Vy) (T 
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(where the cross product is defined by (x x y)'^ = e'^^"' x°^y^). Using that the vector Vx x Vy 
is orthogonal to Vx and Wy, one easily verifies that 



1 



FxFy - - {pxpy + 



VxVy) 



2 



{\PxVy + PyVx\^ - \Vx X Vy\^) . 



This is a polynomial equation for the matrix FxFy, and thus the eigenvalues of the matrix 
are the roots of the polynomial. This proves ()4.13p . 

The identity ()4.14p could be obtained from (j4.13p by a direct calculation using the 
Schwarz inequality and ()4.5p . More directly, we can deduce from (|4.5p that 

detFx = pI- \vxf < , 

and thus 

A+A_ = det {FxFy) = det{Fx) det{Fy) > . | 

This result raises the question whether a fermion system can be reconstructed from the 
local correlation matrices. We first show that, for given px and Vx, one can indeed construct 
a corresponding fermion system. 

Proposition 4.2 Suppose that for given parameters px £ ^ and Vx G the rela- 
tions ^4-4\ \4-5\ ) are satisfied. Then there is a fermionic projector which realizes px and Vx 
as the corresponding local traces and Bloch vectors, respectively. 

Proof. For any x € M, (|4.3p defines a local correlation matrix Fx. According to (j4.5p . Fx 
has at most one positive and at most one negative eigenvalue. Hence diagonalizing Fx by 
a unitary transformation Ux G ^^(2) gives the representation Fx = Ux^DxUx with Dx = 
diag(— a, (3) and a, /? > 0. We define the fermion matrix by 

Ex-^ = \Dx\iUx. 

Then (j4.1ip is satisfied. Furthermore, the relations (j4.3p ensure that the two columns 
of ^' are orthonormal. Hence defining the fermionic projector by (j4.10p gives the desired 
fermion system. ■ 

However, the reconstruction of a fermion system from the local correlation matrices is not 
unique, not even up to gauge transformations, as we now illustrate by a simple example. 

Example 4.3 In the case / = 2 and m = 3, we consider for a real parameter a the family 
of fermion matrices 

As required, the columns of ^ are orthonormal vectors in H. The corresponding local 
traces and Bloch vectors are computed to be 

pi = P2 = 1 , Vi = -V2 = (0, 0, 1)^ , p3 = 0, i/s = . 

These quantities are independent of a. But the fermion systems corresponding to different 
values of a clearly are not gauge equivalent. ♦ 
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Figure 1: Values of the minimal action for different numbers of space-time points 



In view of this example, it does not seem appropriate to restrict attention to the local 
correlation matrices. Instead, we shall describe the fermion system in discrete space-time 
by the fermion matrix (14. Sp . which is determined modulo the transformations (14. 9p . The 
importance of the local correlation matrices lies in the fact that the Bloch vectors will be 
very helpful for illustrating the effect of spontaneous symmetry breaking. 

The variational principle in the critical case can be analyzed numerically as a nonlinear 
optimization problem (for details see Appendix |A]). The minimum of the action as a 
function of the number of space-time points is shown in Figure [1] In the following sections 
we discuss the different systems in detail, with the emphasis on the symmetry structure 
of the minimizers. 



4.2 Two Space-Time Points 



For two particles and two space-time points, the auxiliary variational principle in the 
critical case is solved explicitly in [4l Example 3.2]. To avoid repetitions, we here use a 
different method and work with the local correlation matrices (|4.2p . According to (|4.4p . 
the two matrices Fi and F2 have the form ()4.3p with 

pi+ P2 = 2, vi = -V2 . (4.15) 

Hence the term Vx x Vy in Proposition 14.11 vanishes. 

1 r ^ ^ , ^ 

A± = 7 PxPy + VxVy ± \PxVy + PyVx 



4 L'^ 



In particular, the X± are real. Furthermore, the identity (j4.14p yields that A+ and A. 
have the same sign. Therefore, the Lagrangian (|2.14p simplifies to 



C[A. 



xy\ 



A_ 



PxVy + PyVx 
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Summing over the space-time points and using that 

2 1 ^ 

4 



P2I 



-V1V2 ='■ v^, we obtain 



x=l 



x,y=l 



y=l 



, x=l 



pI 



x,y=l I 



[pI + pI-A 



where in the last step we used the first equation in (j4.15p . Minimizing the action, we 
clearly want to choose v as small as possible. Taking into account the constraint (|4.5p . we 
obtain 

inf5 = max(pi,p2) [Pi + pl - A ■ 

V 

By symmetry, we can assume that pi > P2- Using the first equation in (j4.15p . we get 

inf5 = inf pl (2p? - 4pi + 3) . 

pie[i,oo] 

The minimum is attained at pi = 1. We conclude that for the minimizer the local traces 
and the Bloch vectors have the form 



Pl = P2 
The \± have the values 



1 



and 



1, 



V2 



-Vl 



for ^11 and A22 '■ 
for A12 and A21 ■ 



A+ 
A+ 



1 , A_ 
, A_ 








(4.16) 



(4.17) 
(4.18) 



Hence all points have timelike separation. The degeneracy in M.lSh means that we are 
just in the boundary case between timelike and spacelike separation. 

According to Proposition 14. 2^ there exists a fermion system which realizes (I4.16p . It is 
also easy to verify that the fermion matrix ^ and the corresponding fermionic projector P 
as given by 



/o 

1 







1/ 



p 











V 



(4.19) 



indeed satisfy ()4.16p . However, in view of Example 14.31 the fermion system is in general 
not determined by the local correlation matrices. Therefore, we now need to rely on [U 
Example 3.2], where it is shown that the fermionic projector corresponding to ()4.16p is 
unique up to gauge transformations (|2.15p with U £ Q = U{1,1) x U{1,1). Likewise, the 
fermion matrix in ()4.19p is unique up to the transformations of the form (|4.9p . 

The fermion matrix in (j4.19p shows that the particle corresponding to the first column 
is localized at the first space-time point, whereas the particle of the second column is 
localized at the second space-time point. This configuration is what one would have 
expected in view of the Pauli exclusion principle. To verify that the fermion system ()4.19p 
is permutation symmetric, we choose the outer symmetry group O = S2 = {1, cj}, where a 
is the transposition of the two space-time points. Setting 



U{1) 



U{a) 
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Figure 2: Bloch vectors for three space-time points 



the relations (j2.17p are satisfied. We conclude that the minimizing fermion system for two 
space-time points has maximal outer symmetry, and no spontaneous symmetry breaking 
occurs. 

4.3 Three Space-Time Points 

For three space-time points, the critical action becomes so complicated that we could not 
compute the minimizers analytically. Numerically, we found that the minimizing fermion 
systems are permutation symmetric. Assuming this permutation symmetry, the fermion 
systems can also be treated analytically in closed form. We first discuss the numeri- 
cal results and then give the rigorous analysis of the fermion systems with permutation 
symmetry. 

The minimizers of our variational principle can be computed numerically using methods 
of nonlinear optimization (for details see Appendix [A]). Computing the local correlation 
matrices (j4.21l4.3|) of the minimizers, the resulting Bloch vectors Vx form a plane equilateral 
triangle, see Figure [21 The local traces are all the same and coincide with the lengths of 
the Bloch vectors, 

2 

\vx\ = Px = Vx G {1,2,3}. 

Since any two such equilateral triangles can be mapped onto each other by a rotation in M^, 
such that the numbering of the vertices is respected, the Bloch vectors of the minimizers 
are unique up to the [/(2)-freedom (j4.6l 14. 7() . This suggests that the minimizing fermionic 
projector should be unique up to gauge transformations. Furthermore, any permutation 
of the labels of the vertices of the triangle can be realized by a suitable rotation of R'^, and 
this indicates that the minimizer should be permutation symmetric, i.e. that it should have 
the outer symmetry group O = S3. Indeed, our numerical solution is unique up to gauge 
transformations and is permutation symmetric. We conclude that the minimizing fermion 
system for three space-time points has maximal outer symmetry, and no spontaneous 
symmetry breaking occurs. Furthermore, all pairs of points are timelike separated, and 
thus no non-trivial discrete causal structure appears. 

In order to find closed expressions for the minimizer, we next characterize all fermion 
systems with permutation symmetry using methods of [7j. 
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Lemma 4.4 In the case f 
ing to the fermion matrix 



2 and m = 3, for any i? > the fermion system correspond- 

\ 



( 



1 

7! 



-2sinh'i? 


sinhi? 
"v/S cosh 1? 

sinhi? 
\— -v/Scosh-!? 





— 2coshi9 

-\/3smhf? 
cosh'!? 

\/3 sinh •& 
cosh 'd 



(4.20) 



has the outer symmetry group O = S3. Every fermionic projector with this outer symmetry 
group is gauge equivalent to the fermionic projector corresponding to !i4.20\ ) for some i9 > 0. 



The corresponding local density matrices |^.^ are of the form |^.3D with 

2 
3 



Pi 



P2 



P3 



Vl 




V2 




V3 




where 



- (1 + sinh^ -d) 



(4.21) 
(4.22) 

(4.23) 



Proof. We point out that Definition 12.41 does not imply that the mapping a ^ U {a) is a 
representation of O, and in general this mapping cannot even be arranged to be a group 
representation (this will be illustrated in the proof of Proposition 14. 6p . However, for our 
fermion system and the group 5*3 one can indeed arrange a group representation: In [71 
Proposition 9.4] it is shown that there exists a group representation a 1— > U{(t) of the 
outer symmetry group O = 83 provided that for the set T C O of all transpositions the 
following conditions hold: 

(A) f/(r)2 = 1 for all r G T. 

(B) For all r, r' G T we have the implication 

[T,r'] = l =^ [U{tIU{t')] = 1. 



(C) For all distinct x,y,z S M, 

where Txy denotes the transposition of x with y. 

For the group 5*3, there are no commuting transpositions, and thus (B) is trivially satisfied. 
To satisfy (A) and (C) we proceed as in [71 proof of Theorem 9.1]. We can assume that 
our fermion system is a simple system (see [71 Section 4]). According to [71 Lemma 9.5], 
U{t) is unique up to a phase. By choosing the phase appropriately, we can arrange (A); 
this fixes the operators U (r) up to a sign. Furthermore, we know that (C) holds with a 
general phase, i.e. for all distinct x,y,z £ M, 

^x,y ^y,z ^x,y Ux^z — ^ 1 , 
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where Ux^y = U{Txy). The sign of the phase factor depends on our arbitrary choice 
of the signs of the operators U{t). But up to the sign, the factor e^'^^^^V'^) is uniquely 
defined. Prom the permutation symmetry of our system we conclude that it is a constant 
independent of the space-time points, i.e. 

Ux,y Uy^z Ux,y Ux,z = ±6*'^ 1 for all distinct x, y, z e M . (4.24) 

Multiplying from the right by Ux,z and from the left by Ux,yUy^zUx,y, we get the same 
relation, but with the sign of reversed. We conclude that e^'^ = ±1, proving (C). Hence 
we can arrange that the mapping a U{(t) is a group representation. 

According to [7, Lemma 2.10], H can be decomposed into an orthogonal direct sum 
of definite subspaces, which are irreducible under the group representation a i— > U{a). 
Taking the direct sum of all positive irreducible subspaces gives a subspace C H. 
Introducing /~ similarly, we obtain the decomposition H = /^©/^. Hence these subspaces 
are both three-dimensional, and we can arrange by a gauge transformation that 



Exl+ = r: and Exl~ = I I VxGM. 



o) = (c 

In /- we choose a basis {ex)xGM such that Ex&y = <Ja;y (0, 1)"^ and write a vector 
ti G in components as v = {vx)x£M- Then the right equation in ()2.17p yields that U is 
a permutation up to a phase, i.e. for every v £ I~ , 



'a-{x) 



Using the gauge freedom together with the fact that [/ is a group representation, we can 
arrange that U is of one of the following two form^: 

(i) U{a)vx = v^^x) 

(ii) U{a)vx = sgn{a)v^^x) 

In both cases, the irreducible subspaces of these representations (considered as represen- 
tations on C^, which we can identify via our basis with /~) are 




= C 1 and = {J^y 



On we can proceed exactly as for . Thus our representation U restricted to is 
again of the form (i) or (ii) . We denote the irreducible subspaces by , J2 C . 

Clearly, U acts trivially on Ji := Jj^ © J^", which is an inner product space of sig- 
nature (1,1). Since the image of P must be negative definite, Ji can be occupied by at 
most one particle. On J2 := ® J2 j U acts as two copies of the same 2-dimensional 
irreducible representation. According to Schur's lemma, the range of P| jj must be even. 
Since J2 has signature (2,2), it can be occupied by zero or two particles. Thus the only 
way to build up a fermionic projector of rank two is to occupy Ji by zero and J2 by two 
particles. 



^ An easy method to see that these two representations are not gauge equivalent is to compute the 
determinant of U{g) for an odd permutation a. In cases (i) and (ii) this gives plus and minus one, 
respectively. Clearly, the determinant of U{(t) does not change under gauge transformations U{(j) — > 
VU{a)V-^ with V G t/(3). 
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The representations of U on J2 and J2 are clearly equivalent. But the subtle point is 
that the representation matrices are the same only if in and we choose appropriate 
bases. If U restricted to 1+ and /~ is of the same form (i) or (ii), we can choose the same 
orthonormal basis in J2 and Jg", for example the two vectors 

This basis gives rise to the factorization 

J = C^^y, (4.25) 

where the components of are the coefficients of the basis vectors u and v, and V is an 
inner product space of signature (1, 1). The representation U also factors, 

U\j = Ui„<»l, J = F , (4.26) 

where Uixr is an irreducible matrix representation on C^. According to Schur's lemma, the 
fermionic projector must be trivial on the first factor, 

P\j = lc2®P, 

where p is a projector onto a one-dimensional negative definite subspace of V. Representing 
the unit vector in the image of p as (sinh-i?, cosh??), we obtain the fermion matrix 

^ ^ _ sinh-i? Vx sinh?? 

^ \Ux cosh-i? Vx cosh-!? 

This fermion matrix is gauge equivalent to (j4.20p with = 0. In the case that C/ is of a 
different form on I'^ and I~ (i.e. that U restricted to /"'' is of the form (i) and U restricted 
to /~ is of the form (ii), or vice versa), we can arrange that the representation matrices 
on and /~ are the same by choosing in the basis (n, v) and in I~ the basis (f , —u). 
Then we again have the factorizations ()4.251 14.26P , and we thus obtain the fermion matrix 



Ex'i' 



Vx sinh"!? —Ux sinh"!? 
Ur cosh?? Vt cosh"!? 



This coincides precisely with (j4.20p . 

The local traces and Bloch vectors follow from a straightforward calculation. ■ 

For permutation symmetric systems, we can easily analyze the critical case of the varia- 
tional principle analytically. 

Proposition 4.5 Considering the critical variational principle in the class of fermionic 
projectors with outer symmetry group O = S3, the minimum is attained by the fermionic 
projector corresponding to the fermion matrix ^.20 ) with 1? = 0. The minimizing fermionic 
projector is unique up to gauge transformations. 

Proof. We compute the X± using Proposition 14. II together with Lemma [4.41 For the roots 
of the characteristic polynomial of Axx we obtain 



4 2,4 
-+V'' ±-v 
9 3 



(4.27) 
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Figure 3: Bloch vectors for four space-time points 



For any x 7^ y, the Bloch vectors in ()4.22p have an angle of 120°, and thus 

,2 



A± 



1 



4 



±v 



4 _ 3 
9 ~ 4 



A straightforward computation yields that 

S 



-V +@l V ] V 

3 V 27 M 3 16 



(4.28) 



(4.29) 



According to (j4.23p . we minimize the action under the constraint u > 2/3. A short calcu- 
lation shows that the minimum is attained at v = 2/3, giving the result. ■ 



4.4 Four Space-Time Points: Spontaneous Breaking of the Parity Sym- 
metry 

In the case / = 2 and m = 4, Theorem 12.51 vields that there are no fermion systems with 
outer symmetry group O = S/^. In other words, the permutation symmetry is sponta- 
neously broken. This effect can be seen in our numerics as follows. Similar as for three 
space-time points, the lengths of the Bloch vectors are all equal and coincide with the 
local traces, 

\v^\ = = i VxG {1,...,4}. (4.30) 

The Bloch vectors of our numerical minimizers form a tetrahedron, see Figure [3l By 
labeling its vertices by the corresponding space-time points 1, ... ,4, we give the tetrahe- 
dron an orientation. As a consequence, two such labeled tetrahedra can in general not be 
mapped onto each other by a rotation in M^, as one sees in the two examples of Figure El 
In other words, there are two minimizing fermionic projectors, which distinguish a left 
or right orientation on the tetrahedrons and are not gauge equivalent. Since the parity 
transformation x — > —x reverses the orientation of a tetrahedron, we can also say that 
the symmetry under parity transformations is spontaneously broken by the minimizing 
fermionic projector. The outer symmetry group of a minimizer is merely the alternating 
group A4 of all even permutations, because only those permutations can be realized by 
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rotations in R^. This fact can also be expressed by saying that the minimizing fermionic 
projector spontaneously breaks the permutation symmetry of the space-time points, giving 
rise to an orientation of the space-time points. 

In the next proposition we give closed formulas for A4-symmetric fermion systems 
satisfying ()4.30p . These fermion systems coincide precisely with our numerically found 
minimizers. Using the methods of one could characterize all fermionic projectors with 
outer symmetry group A^. However, since the group A4 and its central extensions have 
many irreducible representations, this analysis would be lengthy. For simplicity, we here 
only give the fermion matrices, but we do not prove uniqueness. 



Proposition 4.6 For every ip 
matrix 



ib27r/3, the fermion system corresponding to the fermion 







with 



^3 



1 



1 
1 

V 1 





V2 



^J2^ 



(4.31) 



has the outer symmetry group A4. The local correlation matrices satisfy \4-30^ - 

The corresponding Bloch vectors are in a tetrahedron configuration. 

The Bloch vectors of the two systems corresponding to if = 27r/3 and ip = —2it/3 have 
opposite orientation. The corresponding fermionic projectors are not gauge equivalent. 

Proof. The group A4 is generated by the two even permutations 



12 3 4 

13 4 2 



and 



12 3 4 
2 14 3 



Hence it suffices to satisfy the relations (|2.17p for these two group elements. Since all 
transformations will leave the first component of E^'^ at each space-time point unchanged. 



we only consider the transformation of the vector ip = (tpi . 
product induced by H is the standard scalar product on C^. 
matrices 



U{a) 



/I 






1 



0/ 



Uir) 



. . . , -04) , on which the inner 
We introduce the two unitary 



(4.32) 





1 







1 




















i 







—i 


0/ 



A direct computation using the trigonometric identities 

I 2Tvi I — 

\^2e^— = ±iV3, 

yields that 



1 



e 3 



^/3eT 



3 



where V are the following C/(2)-matrices: 







V{t) 



V3\V2 -1 



(4.33) 



Substituting (j4.33p into the formula for the fermionic projector ()4.10p . the matrices V drop 
out, and we obtain the left equation in ()2.17p . The right equation is obvious from (I4.32p . 
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We conclude that really is an outer symmetry group of the fermion system with fermion 
matrix (gSH). 

It can be verified by a straightforward computation that the Bloch vectors correspond- 
ing to (j4.3ip form a tetrahedron and satisfy (j4.30p . Alternatively, these facts are verified 
without a computation as follows. Since acts transitively on M, the local traces px all 
coincide. From (|4.4p we conclude that = 1/2. Next, it is obvious from (|4.3ip that the 
matrix has rank one. Using (j4.1ip it follows that the local correlation matrices 
all have rank at most one. This means in the representation (j4.3p that \vx\ = Px- Hence 
the Bloch vectors Vx have non-zero length and are in an ^4-symmetric configuration. The 
only such configuration is a tetrahedron. 

As is obvious from (j4.3ip . the systems corresponding to ip = 2tt/3 and ip = —2it/3 are 
obtained from each other by exchanging the third and fourth space-time point. Since this 
transformation changes the orientation of the tetrahedron, the tetrahedra of the two sys- 
tems cannot be mapped onto each other by a rotation in M.^. However, the local correlation 
matrices of gauge equivalent fermionic projectors are unitarily equivalent, because (|4.2p is 
gauge invariant and depends only on the choice of the basis in P{H). We conclude that 
the fermionic projectors corresponding to ip = ib27r/3 cannot be gauge equivalent. ■ 



The 744-symmetric fermion systems of the previous proposition are also interesting 
from the point of view of representation theory, because they are examples where the outer 
symmetry group has no representation. In order to get a group representation, one needs 
to construct a central extension N of the outer symmetry group. In the next proposition 
we construct this group extension for the systems of Proposition 14.61 This explains why 
in Definition 12.41 it was important not to demand that the mapping a 1— > U{(j) should be a 
representation of O. The following example also illustrates the abstract constructions in [71 
Section 5]. We denote the cyclic group of order two by Z2 (e.g. one can take Z2 = {1, —1} 
with the standard multiplication). By a Z2-extension of ^4 we mean a group N together 
with group homomorphisms i and vr such that the following sequence is exact, 

— > 1.2 ^ M ^ Ai — ^0, (4.34) 

where denotes the trivial group. 

Proposition 4.7 There is no unitary representation U (a) of the outer symmetry group A4, 
which satisfies \2.11^. There exists a ^2- extension M of A4 together with a unitary repre- 



sentation U which describes the outer symmetry in the sense that for all cr G M, 

U{a)PU{ar^=P and U{a)ExUiay^ = E^^^^^^^x) Vx G {1, . . . , 4} , (4.35) 
where tt : M ^ A4 is the projection in ^.34^ . 

Proof. A short computation shows that the conditions ()2.17p determine the matrices 
in ()4.32p uniquely up to a phase. (Using the notion from [7j, this could be expressed by 
saying that the free gauge group is merely the U{1) of global phase transformations.) In 
order to obtain a group representation, we must certainly satisfy the conditions U{a)^ = 1 
and f/(vr)^ = 1. Hence U{a) must coincide with the matrix in (j4.32p up to a sign, 
whereas U{t) could be modified by multiplication with g^^'^*/^. Computing the product 
U{T)U{a), we find that 

{U{T)U{a)f = ±il. 
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Figure 4: Bloch vectors for five space-time points 

If U were a representation, this would be equal to [/((ro")^) = U{\) = 1, a contradiction. 

In order to construct AA, we first note that changing the phase of U changes the phase 
of the matrix V in (j4.33p accordingly. But since V is only a (2 x 2)-matrix, we can fix its 
phase up to a sign by imposing that its determinant should be equal to one. This leads 
us to introduce the matrices 

and to define the matrices U^{t) and U^{a) by the relations 

The four matrices {U^{a), C/^(r)} clearly satisfy the left equation in (j2.17p . They gener- 
ate a group, denoted by M. The injection l in ()4.34p is defined by t : Z2 — > : ±1 1— > ±1. 
Every matrix U £ M clearly satisfies the right relation in (|2.17p for a suitable a £ A4. This 
defines the mapping vr in (j4.34p . Then the relations (14.35P are satisfied by construction. 
Obviously, l is injective, and its image is the kernel of vr. To verify that vr is surjective, one 
represents any group element g £ as a product of powers of a and r. Taking the prod- 
uct of the powers of the corresponding matrices U^{a) and U^{t) gives a matrix U £ N 
with tt{U) = g. ■ 



4.5 Five Space-Time Points: Spontaneous Breaking of the Translation 
Symmetry 

For five space-time points, we obtained the following numerical results. The Bloch vectors 
of the minimizers are unique up to rotations in space and a different numbering of the 
vertices; Figure H] shows a typical example. The lengths of the Bloch vectors coincide with 
the corresponding local traces, 

l^^xl = Px Vx G {1, . . . ,5} . (4.36) 

A new feature, which we found only for five space-time points, is that the local traces, and 
thus also the lengths of the Bloch vectors, have different values. Namely, px ~ 0.3883 for 
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two of the space-time points, whereas px ~ 0.4077 for the other three space-time points. 
The two shorter Bloch vectors he on a straight hne, whereas the other three vectors he 
in a plane perpendicular to this line and form an equilateral triangle. The subgroup of 
rotations which leave the Bloch vectors unchanged is the group S^, corresponding to the 
three rotations about the straight line by the angles and ±120° and the three rotations 
by 180° around the three Bloch vectors in the perpendicular plane. This group is indeed 
the outer symmetry group. Hence the symmetric group ^5 is spontaneously broken to the 
much smaller group O = 5*3. Again, all space-time points have timelike separation, and 
so the emerging discrete causal structure is trivial. 

The spontaneous breaking of the outer symmetry group to the group can be in- 
terpreted as follows. As usual, we call a space-time homogeneous if the outer symmetry 
group acts transitively. Equivalently, one can take the following definition from [3]. 

Definition 4.8 A fermionic projector P is called homogeneous if for any xo,xi £ M 
there is a permutation a : M ^ M with a{xo) = xi and an unitary transformation U 
such that 

P{a{x),a{y)) = UP{x,y)U-^ Vx,yeM. 

We also refer to the permutation cr as a translation in space-time, which translates xq 
to xi. Since the local traces are not all the same, our minimizing fermion system is 
certainly not homogeneous. Thus the homogeneity of discrete space-time is destroyed by 
the fermionic projector. Using a more graphic notion, we also say that the translation 
symmetry of discrete space-time is spontaneously broken by the fermionic projector. 

Using the symmetry structure of the minimizing fermionic projector, we can compute 
the lengths of the Bloch vectors analytically. After a suitable rotation in we can assume 
that the Bloch vectors of the equilateral triangle are proportional to the unit vectors 

ei = , 62 = (^^/2j , 63 

whereas the two Bloch vectors on the straight line are positive multiples of the unit vectors 

6*4 = and 6*5 = 

vJ 

In view of the symmetries of the numerical minimizers, we take the ansatz 

Vx = acx for X = 1, 2, 3 and Vx = (3 Cx for x = 4, 5 

with two parameters a, f3 > 0. Imposing (j4.36p . the relation Tr(P) = 2 leads to the condi- 
tion 3a + 2/3 = 2, leaving us with one free parameter a £ (0, 2/3). Using Proposition 14. H 
we can write the action as a function of a: 

S{a) = 10 - - 18 + 15 - 6 a + 1 . 
8 

A straightforward calculation shows that there is precisely one real minimum at 

a = --[2 + 2Vrij +-(2 + 2Vl7j +-« 0.4077411555 , 

in agreement with our numerical results. In particular, this calculation shows that the mi- 
nimum of the action is indeed smaller than the action of the more symmetric configuration 
with a = p. 
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Figure 5: Bloch vectors for six space-time points, in an octahedral configuration (left) and 
a non-octahedral configuration (right) 




Figure 6: Bloch vectors for eight space-time points, in a cubic configuration (left) and 
other configurations 



4.6 Six to Nine Space-Time Points: Emergence of a Two-Dimensional 
Lattice 

We now discuss our numerical results for more than five space-time points. For the local 
traces and the lengths of the Bloch vectors we always find 

, , 2 
Px = \vx\ = — Vx G M . 
m 

All the space-time points have time-like separation. But the minimizing fermionic projec- 
tor is not permutation symmetric. 

In the case of six space-time points, the most symmetric polyhedron with six vertices 
is the octahedron. The Bloch vectors of length 2/6 which form an octahedron indeed 
have the minimal action S = 2/27, and we found fermion systems which realize these 
octahedral Bloch vectors. However, we also found another type of minimizer which, within 
the accuracy of our numerics of nine digits, has the same value of the action. For this 
minimizer, the angles between the Bloch vectors are a bit different. In Figure [5] the Bloch 
vectors of these two types of minimizers are shown. 

Similarly, for eight space-time points the symmetric polyhedron with eight vertices is 
the cube. The action of this configuration is 1/24, and this is indeed the minimum. But 
the minimizer is not unique. There are several local minima, as is illustrated in Figure [H 
Each minimizer achieves the same minimum up to nine valid digits. 
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Figure 7: Bloch vectors for seven (left) and nine (right) space-time points 



For seven and nine space-time points, the Bloch vectors of typical minimizers are 
shown in Figure [71 These figures reveal the qualitative behavior of two-particle systems 
for many space-time points: For large m, the Bloch vectors of the minimizers have at 
least approximately the same length 2/m and can thus be identified with points on a two- 
dimensional sphere of radius 2/m. The critical variational principle aims at distributing 
these points uniformly on the sphere, i.e. it tries to maximize a certain distance function 
between these points. For many space-time points, the resulting structure looks locally 
like a two-dimensional lattice on a sphere. The qualitative behavior should be very similar 
to the situation described in the survey article [11], where a general class of repulsive forces 
between points on the sphere is considered. 

5 Two-Particle Systems for the Variational Principle with 
Constraint 

5.1 Two Space-Time Points 

As in Section 14.31 we first discuss our numerical results and then construct the minimizers 
analytically. The first question is for what values of k the constraint ()2.2p can be fulfilled. 
Our numerics (see Appendix [B] for details) shows that this is possible for all k > 2. 
For any k, in this range, we found up to gauge transformations a unique minimizer of 
the variational principle (j2.4p . This minimizer has permutation symmetry. If /t = 2, 
the minimizing fermion system coincides precisely with that of the the critical case as 
constructed in Section IT2l In the case k > 2, the systems are different because the length 
of the Bloch vectors is strictly larger than the local trace. All space-time points have 
timelike separation, and thus the discrete causal structure is trivial. 

The variational principle with constraint (12. 4p could be treated analytically. In order 
to keep the analysis reasonably simple, we only consider fermion systems with permutation 
symmetry. Such systems were already constructed in [4]. 

Lemma 5.1 In the case / = 2 and m = 2, for any 'd > 0, the fermion system correspond- 
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ing to the fermion matrix 



/sinhi? \ 
cosh 



(5.1) 



sinh "i? 
\coshT? J 

has the outer symmetry group O = 82- Every fermionic projector with this outer symmetry 
group is gauge equivalent to the fermionic projector corresponding to \5. 1\) for some "i? > 0. 
The local density matrices |^.^ are of the form with 



Pi 



Vl 



P2 



1 



-V2 



1^1 + 2sinh2T9) 



(5.2) 
(5.3) 



-1, 



Proof. In ^ Example 3.2] it is shown that, after a suitable gauge transformation, the 
fermionic projector can be written in the general form 



P 



I — sinh^ a 


sinh /? cosh /3 

\— sinh a. cosh a 





cosh^ (j 



sinh a cosh ol\ 

sinh /3 cosh /? 

- sinh^ /? 

cosh^ a ) 



with a, /3 G M, and by an additional gauge transformation we can arrange that a, /3 > 0. 
For a permutation symmetric system, the local traces p\ and p2 must clearly coincide. 
Computing them, we obtain the condition a = /3. This gives rise to the fermion ma- 
trix (15. ip . The relations (j5.2l I5.3|l follow by a straightforward computation. ■ 



For permutation symmetric systems, we can analyze the variational principle with con- 
straint in closed form. 



Proposition 5.2 Considering the variational principle ^2.4^ in the class of fermionic 
projectors with the outer symmetry group = 82, the constraint can he fulfilled if and 
only if K>2. Then 



min I A 

x,y=l 



2 I 

xy\ 



+ - , 



(5.4) 



and the minimum is attained by the fermionic projector corresponding to the fermion 
matrix i5.1\) with 

1-F2sinh2i? = (k- 1)^ . (5.5) 
The corresponding Lagrange multiplier fi in 112. 6\) has the value 



P 



1 + 



1 



1 



(5.6) 



Proof. Setting v = \vi\ = \v2\ > 1, a straightforward calculation using Lemma l5.ll and 
Proposition 14.11 yields 



for All and A22 ■ 
for A12 and A21 : 



A± = ^(l + ^;'±2t;) = Ul±vf 



4 
A_ 



^{l + v){l-v), 
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and a short computation gives 

^2 \A.y\' = 2{\Au\' + \A,2\') = l + v' 
x,y=l 

E 14,1 = 2{\Al,\ + \Al,\) = \{l + vr. 

x,y=l 

Hence the constraint (12. 2p reduces to the equation 

Hence the condition v > 1 can be fulfiUed only if k > 2, and in this case v = {k — l)^. 
This yields (|5.5p and (j5.4p . To determine the Lagrange multiplier, we first compute the 
action (j2.5p . 

= v^ + {l- fi) (1 + v^) , 
and demand that its u-derivative should vanish at v = {k — 1)^ . ■ 

It is worth noting that the value of the Lagrange multiplier (j5.6p is strictly larger than 
the critical value iJ. = ^. Hence the corresponding action S^, (j2.5l 12. 6p . is unbounded 
below. We thus have an example where a minimizer of the variational principle with 
constraint (j2.4p fails to be a minimizer of the corresponding auxiliary variational princi- 
ple 



5.2 Three Space-Time Points: Emergence of a Discrete Causal Struc- 
ture 

For systems with three space-time points, we can study the emergence of a non-trivial 
discrete structure in detail. We begin by a rigorous analysis of the variational principle 
with constraint for permutation symmetric fermion systems. We will show that for k 
larger than a critical value, space-like separation appears. Afterwards, we will discuss our 
numerical results, which confirm the analytic picture, except that for very large values 
of K in addition the permutation symmetry is spontaneously broken. 

The permutation symmetric fermion systems were characterized in Lemma|131 making 
it possible to analyze the variational principle in detail. 

Proposition 5.3 Considering the variational principle 1^2. 4\ ) in the class of fermionic 
projectors with the outer symmetry group O = S3, the constraint \2.2\) can he fulfilled if 
and only if k > | . Then the local correlation matrices are of the form with 



— and \vr\ 
3 ' ' 



The function 12. 3\) takes the value 



^(72ac-32)3 



ifK < 

- {12 + 6V-32 + 81k)^ ifK > 



68 
81 
68 
81 



Z[P] 



- (V18K-8 + 1) 

y 



if K < 

- 81 

8 , , , K , 68 

-(2 + V-32 + 81k) + - ^fK>- 



(5.7) 
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In the case k < ^ all space-time points {x,y) are timelike separated. If conversely k > H, 
different space-time points are spacelike separated. 

Proof. In the proof of Proposition 14.51 we already computed the roots of the closed chains 
()4.27|4.28p of a fermionic projector with Sa-symmetry. The closed chains A^x always have 
the same spectral weight: 

\Axx\ = 7^ (4 + 9v^) for all x G M. 
18 

Regarding different space-time points x ^ y there are two cases to look at. If w < ^ the 
closed chain A^y has two real roots and therefore 

\Axy\ = ^(8-9^^^) for all x^y. 
The constraint P £ V{k) gives the condition 

32 + 81t;^ = 72k, 

which is fulfilled exactly in the case v = |(72k — 32)4. The condition v < is therefore 
equivalent to k < H. 

In the remaining case v > we obtain a pair of complex conjugate roots and therefore 

1 



1^^3,1 = 2v^A+A_ = —(4 -91-2) for ah x^y. 

io 



The constraint P G 'P{k) gives the condition 
which is fulfilled if and only if 



243^'' - 72v^ + 48 = 108k, 



Writing Z as 



V = -(12 + 6V-32 + 8lK)i 
9 



and using (|4.29p gives ([5 

Comparing the values of Z of the above fermionic projectors to the minimal values cal- 
culated numerically, one finds that the fermionic projector with full Sa-symmetry is the 
minimizer for small k. If however k > 1, we found fermionic projectors numerically which 
realize even smaller values of Z. This is shown in Figure [U where the curve corresponds 
to the analytic minimizers with Ss-symmetry, whereas the single points are the numerical 
minimizers. Since these points lie below the curve, the corresponding fermionic projec- 
tors cannot have Sa-symmetry. Hence for the minimizers for large k, the permutation 
symmetry is spontaneously broken. 

In the case k > H , the discrete causal structure is also a causal set (see for example ) , 
albeit in a trivial way where each point has timelike separation only from itself. We point 
out that the discrete causal structure of Definition 12.61 differs from the definition of a causal 
set in that it does not distinguish between future and past directed separations. Thus in 
order to obtain a causal set, one would have to order the points in each time- like separated 
pair in such a way that the resulting partial ordering is transitive. It is an open question 
whether this can be done for general fermion systems in discrete space-time. 
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Figure 9: Plots of A+ and A_ for varying k 



A Numerical Implementation of the Critical Case 

We first bring the critical Lagrangian (I2.14P into the form which is most convenient for 
the numerical analysis. The roots of the characteristic polynomial of a 2 x 2-matrix Axy 
are 

A± = ^ {TT{A.,y) ± , 

where 

^xy — Tr(j42;y) 4det 

If /S.xy < 0, the X± form a complex conjugate pair, and the Lagrangian vanishes. If 
conversely A^y > 0, A-|- and A_ are both real and according to ()4.14p have the same sign. 
Thus in this case the critical Lagrangian (l2Ti]l simplifies to £ = (A+ - = A^.y/2. 

We conclude that ^ 



31 



where is the Heaviside function, and thus 

x,y£M 

We are now making use of the gauge freedom. For every non-vanishing 2 x 2-matrix there 
exists a unitary U £ U{1,1) which transforms the matrix in exactly one of the following 
forms: 

v\ f u v\ fuwi 

or or 

U W J \U W J \U W2 

where u,v > and w,wi,W2 € C. We now let ^ = (V'i)V'2) be a fermion matrix for two 
particles, so 

<ipi I Vi> = -Sij . (A.l) 

We can assume that after a suitable gauge transformation, every local fermion matrix 
is of the first form: 

v^^ 



Ex^ = { ^ , (A.2) 

\Ux WxJ 

where Ux,Vx > 0, and Wx = Xx + iyx £ C. This is justified as follows. For two space-time 
points, the explicit calculation in Section [4.21 shows that the minimizing fermion matrix is 
of the form ()A.2p . For three space-time points, we studied all possible cases systematically 
and found that the ansatz (1A.2P gives rise to all minimizers. Finally, for more than three 
space-time points we also did numerics with an alternative algorithm which does not use 
the gauge freedom. This algorithm, which is less efficient, always gave the same minimizers 
as the ansatz ()A.2p . 

A direct calculation using (jA.2[) and (j4.10p gives the fermionic projector and the closed 
chains. Then the Lagrangians C[Axy] map the vector 

C — (^1 ; • • • ) ^m; ^1 ; • • • ) 'Vm i ^1 1 • • • ? iVli • • • i Um) 

to a real positive number, and summing over all space-time points gives the action S. S is 
the function to be minimized, where the vector ^ has to lie in the feasible set given by ()A.ip . 
This leads to four constraint functions mapping M^™ to M. Namely, the normalization 
constraints are realized by the functions 

m m 

n(0 = ^? - 1 and r2{0 = Y^^xj + yf - vj) - 1 , 

i=l 1=1 

whereas the orthogonality condition gives rise to the two functions 

m m 

1=1 1=1 

Our optimization problem is to 

minimize S{^) where ri{^) = for i = 1, . . . ,4 . 

We solve this nonlinear optimization problem with constraints with the penalty-method 
(see [To], Chapter 17), where each minimization step is carried out with the algorithm of 
Fletcher- Reeves (see Chapter 5). Although the critical action is only piecewise C^, 
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working with the gradient of S works well, probably because the action is smooth near 
the minimizers. We choose the quadratic penalty-function Q to be 

4 

1=1 

More precisely, our algorithm proceeds as follows. Starting with a random vector ^g, the 
penalty parameter Lq = 1000 and the tolerance tq = 10~^, we minimize in each step 
Q(.;Lk) for a fixed parameter L^. For this minimization we use Fletcher-Reeves, starting 
at and terminate the algorithm if ||V(5(Cfc; < Tk or if the loop is ran too often, 
i.e. 100000 times. If the calculated minimizer is feasible, i.e. if Yli=i '''i i^k) > 10^^'^, 
the algorithm stops. Else the penalty parameter is increased to -L/t+i = 1-1 Li^. and the 
tolerance is reduced to Tk+i = 0.9 Tk- The calculated minimizer S^k+i becomes the starting 
vector in the next step, £,k+i ~ ^k- The C-program leads to the following minimal actions 
depending on the number of space-time points m: 



m 


1 


2 


3 


4 


5 


S 





1 


1/3 


1/6 


0.10701459. . . 


m 


6 


7 


8 


9 


10 


S 


2/27 


0.05442177. . . 


1/24 


0.0329218. . . 


2/75 



B Numerical Implementation of the Variational Principle 
with Constraint 

We now describe the numerical implementation of the variational principle with con- 
straint (12. 4p . A major difference compared to Appendix E] is that we do not use a gradient 
method. Namely, gradient methods did not work, because the function Z does not seem 
to be smooth near the minimizers in the class V{k). 

The smallest k for which the set is non-empty, denoted by Kmin, is computed 

with the following algorithm. 

Step 1: choose all entries of an initial fermion matrix, i.e. the vectors ui, . . . ,Uf, randomly 
in C^*", where the absolute values of all real and imaginary parts are smaller than 
1; introduce an initial distance 6, set S = 1 

Step 2: calculate the fermionic projector P and for all x,y £ M the closed chains Axy 

Step 3: calculate 

with the penalty function 

/ /-I / 

^(■Ui, . . . = Lnorm ^ |<Wi|Uj> + 1| + Lorth ^ ^ \<Ui\Uj>\ , 

i=l i=l j=i+l 

where i^norm and Lorth ai'e positive parameters which can be increased stepwise while 
running the algorithm; for the minimizer of Kmin the value of the penalty function 
is zero (i.e. all normalization conditions are fulfilled), and we get Kmin 



= ^ \Axyf + a{ui,...,Uf) 

x,y£M 
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Figure 10: Minimal value of k for various fermion systems 



Step 4: in order to get neighboring fermion matrices choose increment vectors Aui, . . ., 
Auj, multiply them by 5] the performance of the algorithm increases with the num- 
ber of considered neighboring fermion matrices; this number is to be chosen depend- 
ing on the available CPU power and the size of the considered fermion system. 

Step 5: calculate Hjain for ui + Aui, . . ., u/ + Auj, there are two cases: 

Case 1: If there is a fermion matrix with a smaller Hjaim then consider it as the 
new initial matrix and go back to step 4. 

Case 2: If none of the considered nGigliboring fermion matrices has a, smaller /^min) 
go to step 4 with a smaller (5, for example put 5 



STOP when 5 < 10" 



We considered fermion systems consisting of one to five particles and up to ten space- 
time points. The numerical results are shown in Figure [101 For one- and two-particle 
systems we find that the Kmin coincide precisely with the minima of the critical variational 
principle. Furthermore, for all considered fermion systems we observe that 5min('^min) = 

After having calculated Kmin for a given fermion system, we now solve the variational 
principle for any given k > Kj^m- We use the same algorithm with one exception: instead 
of minimizing Huiin, we minimize 



Z(k) 



x,y£M 



I I 

x,y£M 



+ a{ui, . . . ,Uf), 
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Figure 11: Minimal Z for variations in V-^{k) (x) and in V{k) (+) for the system with 
three particles and four space-time points 



where Lgide is another positive parameter. The additional penalty term 



|2 



|2 



K. 



gives rise to the additional constraint equation YlxyeM \ ^xy\ 

We finally mention some numerical results which are not referred to from the main 
sections. A disadvantage of the above algorithm is that the normalization conditions 
are difficult to fulfill, and thus we had to discard many fermion configurations (for two 
particles and three space-time points, for example, we worked with 10 million neighboring 
matrices). Therefore, it is worth considering whether the normalization conditions can be 
relaxed without changing the minimizers. To study this question, we introduced a new 
variational principle with a more general operator P. This new variational principle is 
solved numerically, and we compare the results to the original variational principle with 
constraint. 

Definition B.l A self-adjoint operator A in an indefinite inner product space is called 
positive if <u\Au> > for all u ^ H. 

For a fermionic projector we have Tr(P) = /, the rank of P is equal to /, and (— P) is 
positive due to the relation 

<u|(-P)-u> = -<u\P'^u> = -<Pu\Pu> > 

Therefore, the next definition extends the class of fermionic projectors. 

Definition B.2 An operator P in an inner product space {H, <.|.>) is of the class pf , 

if {—P) is positive, Tr(P) = f and rk(P) < /. 

For a given parameter k > we define a family of operators 

-pf (k) = Ip of the class P^ with ^ \A^y\^ 

I x,yGM 
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Figure 12: Minimal Z for variations in V^{k) (x) and in 'P(k) (+) for the system with 
four particles and four space-time points 

We now introduce the modified variational principle 

minimize Z = \^xy\ varying P in the class V-^{k), 

x,y&M 

where the number of particles / and the discrete space-time {H^ <•!•>, {Ex)x&m) are kept 
fixed. For this variational principle, the existence of minimizers is proved in Theo- 
rem 2.8]. 

In order to solve this variational principal numerically, the existing algorithm must be 
modified. In step 1 for i G {!,... , /} any vectors Ui with <Ui\ui> < are admissible. 
This ensures the positivity of (— -P). The condition Tr(P) = / is realized by a simple 
rescaling of the vectors Ui. The penalty function is no longer needed, so here we choose 
a{ui,...,Uf) = 0, which reduces the number of numerical operations and makes the 
algorithm faster. 

Again we first calculate Kmin which is the smallest possible k. We find for all consid- 
ered fermion systems that the value of Kmin coincides with that for the original variation 
principle. In Figure [11] and Figure [12] we compare the minimal values of Z for operators 
in the class (points x in green) with the corresponding values for projectors (points + 
in red). In general, the minimum in the class is strictly smaller, which was expected 
because we minimize over a larger class of operators. However, it is remarkable that 
for K = Kmin the minima are the same. It turns out that the minimizers in the class V-^ 
are actually projectors. Hence in this case, the minimizers in automatically satisfy the 
normalization conditions. 

Acknowledgments: We thank Joel Smoller for helpful comments on the manuscript. 
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